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Abstract. Diluted magnetic semiconductors (DMS) like Gai-^Mn^As are described by a realistic tight- 
binding model (TBM) for the (valence) bands of GaAs, by a Zener (J-)term modeling the coupling of the 
localized Mn-spins to the spins of the valence band electrons, and by an additional potential scattering (V-) 
term due to the Mn-impurities. We calculate the effective (Heisenberg) exchange interaction between two 
Mn-moments mediated by the valence electrons. The influence of the number of bands taken into account 
(6-band or 8-band TBM) and of the potential (impurity) scattering V-term is investigated. We find that 
for realistic values of the parameters the indirect exchange integrals show a long-range, oscillating (RKKY- 
like) behavior, if the V-term is neglected, probably leading to spin-glass behavior rather than magnetic 
order. But by including a V-term of a realistic magnitude the exchange couplings become short ranged 
and mainly positive allowing for the possibility of ferromagnetic order. Our results are in good agreement 
with available results of ab-initio treatments. 



1 Introduction 

Dilute magnetic semiconductors (DMS) have risen a lot 
of interest during the last decade because of their possible 
use for spintronic devices [JJ (e.g. as spin field effect tran- 
sistor or for spin-polarized light-emitting diodes). A key 
goal is the achievement of Curie temperatures Tc well 
above room temperature; the largest Tc obtained so far 
is, however, only of the magnitude of about 173 K and has 
been obtained for Gai_^Mn^As [2J. The possibility of a 
much larger Tc has been predicted theoretically, and even 
a T c « 700 K has been predicted for GaMnN 0, but this 
could not yet be confirmed experimentally. Furthermore, 
it has been demonstrated that the theoretical prediction 
of such a huge Tc is the consequence of several drastic 
approximations, namely (1) perturbative treatment, (2) 
mean-field approximation and (3) virtual crystal approxi- 
mation for the treatment of disorder |4l5j . Therefore, im- 
proved and more reliable theories of DMS are necessary, 
and so far two different kinds of approaches have been 
applied, namely (1) realistic band structure model stud- 
ies (kp-, Kohn-Luttinger, empirical tight-binding models) 
[5] and (2) material specific ab-initio calculations [7]. Re- 
cently it has been shown that the essential properties of 
DMS can already be described within simplified model 
studies [8 9 using a one-band model. Therefore, a model 
with a more realistic band structure and a non-perturbative 
treatment of the exchange coupling between the carriers 



and the magnetic moments should be able to realistically 
capture already the essential physical properties of DMS. 

In this paper we start from a realistic empirical tight- 
binding model (ETBM) for the electronic properties of 
GaAs and take into account the Mn-impuritics by two 
additional terms, namely the local exchange coupling J p( [ 
between the spins of the (valence p-) electrons (or the 
holes) and the localized magnetic Mn-(d-shell) moments 
(J-term) and the impurity potential (V-term) provided 
by the Mn-ions. A similar model was also used in recent 
Monte-Carlo (MC) simulations [TIT]. However, in this pre- 
vious work only the three valence bands (per spin direc- 
tion) of GaAs were taken into account, furthermore the V- 
term was completely neglected and the whole investigation 
was restricted to relatively small systems. In the present 
paper we also start from a realistic ETBM description 
of the electronic properties of GaAs. The TB-paramctcrs 
are determined so that the most important known band 
structure features (energy gap, effective masses, Luttinger 
parameters, etc.) are well reproduced. We consider both, a 
six-band model (i.e. only the three valence bands per spin 
direction as investigated in [TO] ) and a more realistic eight- 
band model (including also one conduction band) . We take 
into account a finite concentration x of Mn-ions, which 
provide for a local impurity potential V, a local exchange 
coupling J pc i and simultaneously serve as dopants (ac- 
ceptors), and treat the disorder exactly. Here the kernel- 
polynomial method (KPM) [IT] is employed, which allows 
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for the essentially exact calculation of the density of states 
of the disordered system for relatively large (finite) sys- 
tems. The KPM can also be used for the calculation of 
the effective (long-ranged, Heisenberg) exchange interac- 
tion Jij between two Mn-ions at sites i and j mediated by 
the free carriers (i.e. the electrons or rather the holes in the 
valence bands), starting from the Lichtenstein-Katsnelson 
formula |12j . 

For realistic choices of the parameters {J p dS = 3 eV, 
Mn-concentration x = 5%) we obtain a RKKY-like, long- 
ranged oscillating behavior of Jij, if the potential scat- 
tering V-term is neglected, for both, the 6-band and the 
8-band model. In particular, Jy may be positive or neg- 
ative (ferromagnetic or antifcrromagnetic) depending on 
the distance between two Mn-ions and on the configura- 
tion. Therefore, no magnetic order can be expected but 
rather a spin-glass behavior. On the other hand, if the V- 
term is properly included (of the correct magnitude so that 
the bound state due to Mn impurities in the gap is repro- 
duced) we obtain a more short ranged Jy , which is mainly 
positive so that ferromagnetic order should become possi- 
ble. Altogether the 8-band description is, of course, more 
reliable and realistic and leads, in particular, to valence 
bands of a realistic band width in the simplest ETBM de- 
scription. We also demonstrate that for the 8-band model 
and a finite V our result for Jij is already in rather good 
agreement with the corresponding results of ab-initio cal- 
culations. 



2 Theory 



Our approach to theoretically describe the physics of the 
DMS material Gai-^Mn^As combines an empirical tight- 
binding model (ETBM) for the host semiconductor with 
additional terms describing the nonmagnetic impurity scat- 
tering and a spin-spin interaction term, which couples im- 
purity spins to carrier spins. This kind of model is referred 
to as the V-J model in the literature, see Refs. |13ll"4"] . but 
we want to go beyond a single-band description. Neglect- 
ing the potential scattering V our ansatz becomes equiva- 
lent to the Zener model. From the corresponding Green's 
functions one is able to calculate the effective Mn-Mn ex- 
change integrals in the presence of disorder according to 
the Lichtenstein-Katsnelson formula, Ref.|12]. We treat 
the disorder exactly (for finite systems) and perform the 
calculations for different realizations of the disordered sys- 
tem and, therefore, obtain a distribution of and calculate 
the average of the couplings Jij . These couplings can then 
be used as the exchange interaction of a disordered Heisen- 
berg model to calculate magnetic properties of the system. 
After this short overview about the methodology applied, 
we are going to be more specific. 



2.1 Multiband V-J model 

We study the following hamiltonian: 

(1) 

Here, the first term corresponds to the ETBM description 
of the relevant energy bands of GaAs with tf? being the 
hopping integral of an electron from lattice site j to i in 
the orbital /3 to a. The corresponding creation (annihila- 
tion) operator is c| (cjpa) with the additional spin- index 
a 6 {1\4}- At this point it is convenient to discuss the 
tight-binding parametrization for the t"? and the basis in 
more detail. One possibility is to set up the hamiltonian 
for the full atomistic lattice and use s,p, d-orbitals be- 
cause the chemical bond of GaAs is sp 3 -hybridized and the 
Mn 2+ impurities under consideration introduce d-orbitals 
into the problem. These kinds of approaches were already 
studied in the work of Tang and Flatte [TS] without d- 
orbitals, by Masek et al. 16 including the <i-orbitals as 
well compared in Ref. [17] by Turek et al. We use an al- 
ternative, simpler version of an ETBM starting from well 
localized (Wannier) orbitals on the fee lattice and can use 
two different analytical parametrizations for the t^ . The 
first one uses only p-orbitals describing the valence bands 
of GaAs and was already studied in Ref. (Til], while the 
second one uses the sp 3 -basis as first suggested by Loehr 
in Ref.[TB] and was also applied to ferromagnetic semi- 
conductor superlattices [Tj5]. The ETBMs which operate 
on a Bravais lattice are usually referred to as effective- 
bond orbital model (EBOM) and have the benefit, that 
they do analytically include the k ■ p-hamiltonian in the 
limit k — > 0. This means, in particular, that the t"? are 
given by analytic expressions in terms of (anisotropic) ef- 
fective masses m c , Luttinger parameters ji, bandgaps E g , 
critical point energies, the Kane parameter E p , etc. We 
use the input parameters given in Tab.([TJ for the 6-band 
and 8-band ETBM, but neglect the spin-orbit splitting. 
There is no a-priori justification for doing this, since for 
GaAs A so = 341 meV is relatively large. Nevertheless we 
want to compare our results for the effective exchange in- 
tegrals Jij to results from ab-initio tight-binding linear 
muffin-tin orbital theory (TB-LMTO) with the coherent 
potential approximation (CPA) which also neglect spin- 
orbit interaction, see Ref. |20) for details. 

The second term in Eq.([T]) describes the nonmagnetic 
impurity scattering with the short-range on-site potential 
V a and fiiaa- is the occupation operator. For the simula- 
tion of the disorder the variable pi takes the value 1 with 
probability x and with probability 1 — x where x de- 
notes the impurity concentration of the sample for each 
lattice site i. This contribution to the hamiltonian does 
also arise from the Schricffer-Wolf transformation [22] of 
the Anderson hamiltonian, Ref. [33] . We use V a — V for all 
orbitals from now on in order to keep the model param- 
eters at bare minimum and its magnitude will be physi- 
cally motivated and discussed in the results. Of course this 
part of the hamiltonian could be augmented by a nearest- 
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Table 1. Model parameters for GaAs used in the calcu- 
lations. Unchecked quantities (-) are not included in the 
respective analytical tight-binding parametrization. 



Parameter 


Value 




6-bands 


8-ban 


a 


A 


5.64 


032 




V 




eV 


0.341 


£□ 






m c 


m 


0.067 


HP 




V 


-Tic - r±5 V 


eV 


1.519 


£H 




V 




eV 


2.1 






V 


Xsv 


eV 


-3 






V 


X3 V 


eV 


-7 






V 


7i 




6.85 


032 




V 


72 




2.1 






V 


73 




2.9 


m 




V 



neighbor Coulomb potential similar as done in Ref . [24j . 
but we choose to study the on-site case solely because it 
prevents introducing even more parameters in our study. 

The last term of Eq.([T]) couples impurity spins Si to 
carrier spins Si jQ of the host system via a local contact- 
interaction with coupling constant J a and we have to deal 
with quantum mechanical spin 5/2 operators for the Mn 2+ 
impurities and 1/2 for the carriers (holes). In the case of 
Gai-zMn^As the fundamental magnetic coupling mecha- 
nism is the coupling J p d between the spins of the (p-like) 
valence electrons (holes) and the impurity spins (formed 
by the Mn-d-shells) ; of course, there is also a coupling 
J s d between the (s-like) conduction electron spins and 
the Mn-spins. The magnitude is about J p d » 1.2 eV and 
J s d ~ 0.02 eV, see Rcfs.|25 26 , and we decided to use 
J a = J P d for all orbitals which is again due to our con- 
straint to keep the model paramaters at minimum. Now we 
must explain how we treat this spin-spin interaction term 
in our approach. Our calculations are performed for T = 
K and thus transverse impurity spin-fluctuations should 
be absent. We must then only consider the z-component 
of the spin-operators and replace Sf — > (§*} = S by its ex- 
pectation value for the calculations of the electronic prop- 
erties only. The spin-spin interaction term is then given 

M& = 1): 

^-^2Pi(ni,aA-n it a,±). (2) 

Finally we are left with a single-particle hamiltonian de- 
scribing non-interacting electrons on a lattice via kinetic 
energy, magnetic and nonmagnetic scattering which we 
have to solve for each random configuration of substitu- 
tional disorder. 

In order to get an impression of the quality of the 
parametrization for bulk GaAs, we have plotted the band- 
structure along the corresponding irreducible wedge for 
the 8-band model in Fig.((T]). Due to the analytical expres- 
sion for the hopping matrix elements, the experimental 
bandgap and effective masses at the -T-point are repro- 
duced, see Ref. [21] for an overview, as well as the critical 
point energies at the A-point can be adjusted. In compar- 
ison to first-principle bandstructures in e.g. Refs.[2?28 



Bandstructure of GaAs 




Irreducible wedge 



Fig. 1. Bandstructure of bulk GaAs calculated within the 
8-band model without spin-orbit splitting. 



the ~12 eV bandwidth of the 8-band model is in reason- 
able agreement and the position of the A-point energies 
lies within the values reported in the literature, see e.g. 
Ref. 27 for a comparison. Deviations in terms of energetic 
position and curvature start to appear for larger fc-values 
along the W—L path and along X — W — K as expected. In 
the interval L — T — X the agreement is better, despite the 
first-principles calculation in Ref. [37] has difficulties to re- 
produce the experimental bandgap. The bandwidth of the 
6-band model is approximately «30 eV and the agreement 
is only reasonable in the vicinity of the -T-point, since the 
A-point energies are not pinned. 



2.2 Limitations of the model 

The TB-parameters used in our calculations are chosen so 
that the pure GaAs band structure is reasonably repro- 
duced and that the Mn-impurities are properly described. 
For the latter purpose the impurity parameters are chosen 
so that the bound Mn impurity state in the band gap of 
GaAs is reproduced, and these parameters (chosen for the 
single Mn impurity within GaAs) are used also for finite 
concentrations x of Mn. One should be aware of the fact 
that this treatment is justified only in the low concentra- 
tion limit of small x. Though in principle the method can 
be applied to arbitrary x (even to the case x = 1 of pure 
MnAs, see e.g. Ref. [29]). the parameters are not appropri- 
ate for too large x. One has, for instance, to choose and 
fit also off-diagonal hopping matrix elements appropriate 
for MnAs if interested in higher Mn-concentrations x 



2.3 Evaluation of exchange integrals 

Having solved the electronic problem, one is able to cal- 
culate the effective Mn-Mn exchange interaction between 
two impurities located at different lattice sites and Yj 
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according to the expressions derived in Ref.(12|: 

oo 

J y = i^2 5 / /(^)Tr Q {4(^)^(^)4.(^)^( w )} dw - 

— OO 

(3) 

Thereby vertex corrections have been neglected, which is 
justified according to Ref . [50] . Here, f(u) is the Fermi- 
function, Tr a denotes the trace over orbital indices and 
G^.(oj) is the one-particle Green's function matrix with 
respect to the orbital basis in the spin-sector a. Accord- 
ing to the approximations outlined earlier, the self-energy 
reduces to S . — 3 p dS-l, so that Eq.Q can be written in 
a form where only products of two Green's functions have 
to be evaluated: 

oo 

J « = E ^ 3 / /MGjf V)Gjf (4) 




Via the Fermi function Jy depends on the position of the 
Fermi energy, which has to be determined for the given 
number of carriers, which again depends on the number 
of holes rih introduced by the Mn 2+ -doping. Though, in 
principle, the hole concentration could be less than the 
Mn-concentration x (because of anti-site effects, i.e. Mn- 
ions occupying As-sites), we will assume = 1 for each 
Mn-impurity, because otherwise also additional anti-site 
impurity scattering terms would have to be included. 

The Jij according to Eq.flUl which depend on the ac- 
tual configuration of the disordered system and are, there- 
fore, subject to a probability distribution function, can 
be used as input parameters for a disordered Heisenberg 
model: 

Hueis = — /IPiPjJijS&j. (5) 

i¥=3 



This integral can then be calculated within the KPM for 
all impurity indices i,j and orbital indices a, j3. To give 
the non-expert the basic idea of the KPM we briefly sum- 
marize the concept. The Hamiltonian H is approximated 
in an infinite series of Chebyshev polynomials T n {H) of 
order n defined on the interval [—1,1] and thus the spec- 
trum of H has to be linearly remapped. In practice this 
series must be truncated at some point and the respec- 
tive recursion relation (not shown) is evaluated for the 
moments fj, n (expansion coefficents) up to a certain or- 
der. By using a cosine discretization for the energy-axis, 
the expanded function can be reconstructed by a fast 
Fourier transform of the fi n and in order to prevent so 
called Gibbs-oscillations of the expanded functions it is 
necessary to include a multiplicative kernel, too. We use 
the Jackson kernel as suggested in Ref. II] resulting in a 
Gaussian representation of peaks in the spectral density. 
The broadening S n (ui) of these peaks is analytically con- 
nected to the order of expansion and in addition energy 
dependent, which results in an unphysical inhomogeneous 
broadening not connected to any quasiparticle lifetimes. 
A reasonable choice for the correct order n of expansion 
might be given by the condition, that the average broad- 
ening {8 n (ui)) is smaller than half the energy difference 5e 
of the spectrum of H . Of course one can choose n less or 
higher depending on the purpose, e.g. producing a smooth 
density of states or resolving every spectral peak by itself. 

In our calculations for the Heisenberg exchange inte- 
grals we have checked (not shown) the needed order of 
expansion for a bulk GaAs system with N x = N y = N z = 
N = 16 unit cells in each spatial dimension and decided to 
use approximatively half the order recommended because 
the numerical calculations are still very demanding on a 
cluster using a highly parallelized implementation. In dis- 
ordered systems we would expect an even higher number 
than recommended. For calculations concerning the den- 
sity of states (DOS) we chose the number of moments fi n 
so that a smooth function results. 



2.4 Kernel polynomial method 

We want to study the electronic problem rigorously (for 
a finite system), thereby treating the disorder (impurity 
scattering) exactly. Since the exact diagonalization scales 
with the third power of the dimension of the Hamiltonian, 
we employ the kernel polynomial method (KPM), which 
essentially scales linear with the dimension of the Hamil- 
tonian. The KPM, see the review by Weisse et al. [TT] . 
can be used to calculate the total density of states of the 
disordered system. It can also be used to calculate traces 
over products of matrix elements of operators. Therefore, 
also the Heisenberg exchange interaction in Eq.© can be 
calculated according to: 

oo oo 

[*» [toi^m. 

3 4 J J oj-oj' 

— OO —oo 

t \S(u - H)\ia t) • (ia I \8(u>' - H)\jj3 I)}. (6) 



3 Results 

3.1 Electronic properties 

First we present results for the disorder averaged, spin- 
resolved total density of states (DOS) of the valence bands 
(VBs) for Gai_ x Mn x As with x=5% with and without the 
nonmagnetic on-site scattering term V in Fig|5] The mag- 
nitude of V is chosen to correctly reproduce the acceptor 
level in the gap for Mn 2+ of E b w 113 meV, Ref. 31 , for 
one single impurity in the GaAs matrix which is achieved 
by setting V = 1.93 eV for the 6-band and V = 0.85 eV 
for the 8-band ETBM. The numerical calculations were 
performed for a system of N = 40 conventional unit cells 
in each direction (i.e. 4-N 3 =256 000 lattice sites altogether 
because of 4 sites per conventional unit cell in the fee lat- 
tice) with 2048 moments and the disorder average was 
stopped after the fluctuations were less than 1%. 

Let us start by discussing the 6-band case with V = 
eV; one can observe, that there is no preformed impurity 
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V=OeV V=1.93eV 




-3 -2-1 1 -3 -2-1 1 



V=0 eV V=0.85 eV 




E-E f (eV) E-E f (eV) 

Fig. 2. Disorder averaged density of states (DOS) of the valence bands calculated for Gai-^Mn^As with x—5% within 
different tight-binding models with and without the nonmagnetic on-site scattering term V. The solid lines correspond 
to the total DOS while the spin-up (black) / down (red) sector is plotted in the upper and lower part respectively. The 
system size is N=40 (256 000 lattice sites) and the position of the Fermi-level is denoted with the blue dashed-dotted 
vertical line. 



band (IB) present but a nonvanishing DOS for E — Ef > 
in combination with a small spin-splitting. Furthermore, 
the Fermi-level Ef lies within the valence band. In con- 
trast, choosing a finite V = 1.93 eV changes the situation 
completely. Now an impurity band formation can clearly 
be observed around the acceptor level and for x — 5% it 
is merged with the VB. Due to the fact that we do not 
apply any approximations in terms of disorder, there is a 
pronounced sidepeak present which can be attributed to 
the formation of impurity clusters. Please note, that these 
peaks vanish if one uses the coherent-potential approxima- 
tion for disorder. The Fermi-level lies within the IB near 
the acceptor level, but in the spin-down bandgap, and the 
spin-splitting is largely enhanced. From these facts one 
can conclude, that the system could show an insulating 
character if the eigenstates in the impurity band are lo- 
calized. In the 8-band case the situation is qualitatively 
similar, but as we include more bands, the structure and 
the bandwidth of the DOS is much more realistical. 

So far we have analyzed the VBs only and as a next 
step the conduction band is adressed and visualized in 
Fig-© f° r the 8-band model. Here one observes that in the 
case of V — eV the spin-down DOS is enhanced around 
about E — Ef « 2 eV and for V = 0.85 eV this effect 
does not occur. In both cases there is also a pronounced 



impurity band with sidepeaks in the spin-up sector. This 
is a consequence of using J a — J p d and V a = V for all 
orbital indices a, also for the conduction band. In the case 
of V = 0.85 eV the CB impurity band is slightly shifted 
by « 0.5 eV to higher energies as can be seen from Fig[2] 
Thus the former spin-down DOS enhancement for V = 
eV vanishes in the case of V = 0.85 eV and merges into 
the GaAs CB. 

3.2 Effective exchange integrals 

The effective exchange integrals are plotted in Fig|4] for 
the 6-band model and in FigJS] for the 8-band model. 
The numerical calculations were performed for a system of 
N = 16 (16 384 lattice sites) with 4096 moments and the 
disorder average was carried out over 144 different config- 
urations of disorder while each data point is an average 
over 1008 values. Please note, these calculations are seri- 
ously more demanding than calculating a DOS in terms of 
computing time and are reasonably performed on a clus- 
ter using a parallel implementation. We remark, that the 
positional disorder average must be carried out over equiv- 
alent difference vectors = r\, — with respect to the 
symmetry of the underlying lattice. For instance, an incor- 
rect average over the absolute distance R = ||R»j|| would 
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V=0 eV 



V=0.85 eV 



> 



LU 
I 

LU 



0.05 





°- -0.05 



E-E f (eV) 

Fig. 3. Disorder averaged density of states (DOS) of the conduction band calculated for Gai-^Mn^As with x=5% 
within the 8-band model with and without the nonmagnetic on-site scattering term V. The solid lines correspond to 
the total DOS while the spin-up (black) / down (red) sector is plotted in the upper and lower part respectively. The 
system size is N=40 conventional unit cells in each direction (256 000 lattice sites). 



correspond to an incorrect spherical symmetry instead of 
the correct fcc-lattice structure. Furthermore we close the 
system with periodic boundary conditions and thus the 
correct average has to take this into account as well. So 
we use the shortest distance between two Mn-impurities 
at Yi and Vj and difference vectors whose components are 
larger in magnitude than half the numerical box-size N/2 
are periodically remapped. Also, we calculate the Heisen- 
berg exchange integrals only up to this absolute magni- 
tude of R = ||Rij|| < N/2 because equivalent difference 
vectors with R > N/2 can never be realized for a given box 
size N and the average is then affected too and considered 
as not reliable. 

Comparing the results obtained for the 6-band model 
one can observe that with and without the potential scat- 
tering term the exchange integrals are quite different. For 
V = 1.93 eV the couplings are mainly ferromagnetic and 
at distances R < 4a their magnitude is enhanced while 
at R > 4a a damping seems to occur as can be identified 
in the lower plot in FigJ3] where the couplings have been 
rescaled by the RKKY factor R 3 . Contrary, if the on-site 
scattering term is neglected, the effective exchange inte- 
grals are in general lower in magnitude and the envelope 
part seems to oscillate around zero. The long-range nature 
is only visible in the rescaled exchange integrals. Coming 
back to the former discussion, the results for the couplings 
Jij agree well with the previous electronic picture given. 
As the Fermi-level lies in the VB for V — eV and it is 
known from textbook RKKY theory for metals, that one 
should expect a long-range oscillating behavior of the ef- 
fective exchange integrals in this case, our results are not 
surprising. If the Fermi level falls into the valence band in 
a region with delocalized states, RKKY like oscillations 
can be expected. But if the Fermi level falls into an impu- 
rity band, in which the eigenstates may have a more local- 
ized character, quite a different behavior may result. This 
directly reflects itself in the damped more short-ranged ef- 
fective exchange integrals. Now we shall directly compare 



to the results of the 8-band model in Fig.© and we see, 
that the qualitative features are identical. For V — 0.85 
eV the couplings are slightly reduced in magnitude at dis- 
tances R < 4a in comparison to the 6-band case with finite 

V due to the increased weight of the impurity DOS at the 
Fermi-level. In addition, for V = eV the long-range os- 
cillating tail seems to be more pronounced as the envelope 
part is of larger magnitude at R > 4a compared to the 6- 
band model. Overall the picture of the exchange integrals 
is consistent within the 6- and 8-band model. 

At last we want to compare our results to the exchange 
integrals obtained by ab-initio approaches, see Ref.[2TJ] 
for more details, in Fig. ([6]). Here a direct comparison be- 
tween our 8-band model to results obtained within the 
TB-LMTO+CPA theory is made. One observes, that the 
overall shape and structure of the exchange integrals is 
very similiar for distances up to w 3.25a and only the 
magnitude is different as can be seen from the lower plot 
in Fig.®. In the case of V = eV the couplings are 
mainly lower in value and due to the oscillating nature 
for larger distances, the discrepancy becomes visible at 
about 3 — 3.5a. For V = 0.85 eV the agreement seems 
to be better, since the exchange integrals are larger and 
thus more close to the TB-LMTO+CPA results. In par- 
ticular the 2nd, 4th, 5th and even higher shells do almost 
agree in absolute numbers. From the fact that long-range 
RKKY- like oscillations are absent, the 8-band model with 

V = 0.85 eV seems to be a reasonable model to calculate 
effective Heisenberg exchange integrals and to study mag- 
netic properties. Though the 6-band model with V — 1.93 
eV does also have in principle the same characteristics, it 
provides short-ranged ferromagnetic couplings which are 
partially double the magnitude compared to the 8-band 
approach. Thus from a model point of view, the 6-band 
V-J model is quantitatively less accurate. But the gen- 
eral agreement between our model calculations and the 
ab-initio theory is surprisingly good both in a qualitative 
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and a quantitative manner though we applied a lot of ap- 
proximations on the initial hamiltonian. 

3.3 Critical temperatures 

As mentioned, the Jij calculated and presented in Figs. 
4-6 can be used as the parameters of an effective dis- 
ordered Heisenberg model, for which the Curie temper- 
ature Tc can be calculated. The most simple estimate for 
Tc is obtained from the standard mean-held approxima- 
tion (MFA) for classical spins as applied in Ref . |20] . For 
our Jij obtained for the 8-band model and x = 5% Mn- 
concentration this estimate yields values for Tq IF of the 
magnitude between 160 K for V — and 390 K for V ^ 0, 
which is even larger than that obtained for the ab-initio 
Jij from Ref. [20] ( « 290 K for x = 0.05). However, it 
is clear that the MFA overestimates the true Tc and that 
more sophisticated treatments [713*21 1 11317] of the disor- 
dered Heisenberg model have to be applied. From prelimi- 
nary studies in Ref. [33] (using only the 6-band model and 
a smaller system treated by exact diagonalization) we ex- 
pect realistic T^-values to be of the magnitude Tc < 200 
K. 



4 Summary 

In this work the effective Heisenberg exchange integrals 
Jij between Mn 2+ impurities were numerically calculated 
for Gai-zMn-cAs with x—5% by using two different tight- 
binding models for the host semiconductor. A comparison 
on the relevance of a nonmagnetic scattering term V was 



made with the result, that within the 6-band model for 
the valence bands as well as in the 8-band model includ- 
ing the conduction band a neglection of V gives rise to 
a long-range oscillating RKKY-like tail in the couplings 
Jij. In contrast, the inclusion of a finite nonmagnetic on- 
site scattering term reproduces the Mn-acceptor level in 
the limit of one impurity and changes the nature of the 
couplings completely to a short-ranged and damped char- 
acter and yields mainly positive values. This behavior was 
explained by analyzing the total spin-resolved density of 
states, because by choosing V = eV the Fermi-level lies 
in the valence band and for V > eV an impurity band 
formation around the Mn-acceptor level occurred in com- 
bination with the Fermi-level lying therein. In comparison 
to available ab-initio results, our calculated couplings are 
of qualitatively of the same structure and do even quanti- 
tatively agree. Thus it is important for ferromagnetism to 
become possible to include the potential scattering term. 
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